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Abstract 

The decay of the i^~'"i^'~hadronic atom kaonium is investigated non-perturbatively using 
meson-meson interaction amplitudes taken from leading order chiral perturbation theory in an 
approach adapted from that proposed by Oiler and Oset [IS]. The Kudryavtsev-Popov eigenvalue 
equation is solved numerically for the energy shift and decay width due to strong interactions in 
the Is state. These calculations introduce a cutoff ~ 1.4 GeV in 0(4) momentum space that is 
necessary to regulate divergent loop contributions to the meson-meson scattering amplitudes in the 
strong-interaction sector. One finds lifetimes of (2.2ib0.9) x 10~^'^s for the ground state of kaonium. 
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1. Introduction 



As is the case for the hadronic atoms pionium [l| and kaonic hydrogen (2[ the energy shifts 
and decay widths in the coulomb spectrum of the hadronic atom K~^K~, or kaonium j^, 
also depend primarily on the scattering length of the interacting particles, in this case 
and K~ . This scattering length, in turn, is determined by the properties of the KK strong 
interactions and the associated models for the possible structures of the /o(980) and an(980) 
scalar mesons. These proposals include a kaon-antikaon molecular bound state |in|, a 



q q state [12[, or a qq state |13l-|17 



The molecular state option has been explored in some detail in (g], [8-11||. These calcu- 
lations all employ a SUvi'i) x S'?7yi(3) symmetric model Lagrangian density \d\ to generate 
meson-meson interaction vertices via vector meson exchange. A rather different route was 
followed by Oiler and Oset [l8[. These authors take the meson-meson interaction vertices 
from leading order chiral perturbation theory (xPT) as being an appropriate theoretical re- 
alization of low energy QCD [l9|, and solve the resulting Lippmann-Schwinger equation for 
the corresponding T-matrix elements to provide a non-perturbative description of meson- 
meson scattering and reactions. This approach introduces one additional parameter insofar 
as the meson loops occurring in the resulting equations need to be regularized with a high- 
momentum cutoff. Using a three-momentum cutoff of ~ 1 GeV they find good agreement 
with the measured vrvr — )■ vrvr phase shifts and vrvr — > KK inelasticities up to center-of-mass 
(cm.) collision energies of order ^/s ~ 1.2 GeV. The masses and decay widths of the /o(980) 
and ao(980) scalar mesons are also satisfactorily reproduced as complex poles of the relevant 
T matrices. 

In the following we adapt the OUer-Oset [l8[ approach to study the energy shifts and 
decay widths introduced by the strong interactions into the spectrum of kaonium. 



2. Non perturbative scattering amplitudes 

The calculations reported in [18] start from the xPT leading order interaction Lagrangian 



density [18|, l20|, |21| £2 for the pseudoscalar meson octet at low energies to generate tree- 
level chirally symmetric meson-meson interaction amplitudes, or 4-point vertices, of good 
isospin / for s-wave meson scattering. The explicit form of £2 is given by Eq. ( !T7|) of the 
Appendix. These vertices form a symmetric matrix V/j, where the labels i and j each refer 
to the channels {KK}, {tttt} or {tt^t]}. (We retain the labelling convention of [l8| by setting 
ihj) = (1) 2) where 1 refers to {KK} for both / = and 1, while 2 refers to {vrvr} for 1 = 
or {vr°?7} for 1 = 1). 

Knowing the V/j one can construct the coupled integral equations for the associated 
scattering matrices T^y However, for s-wave scattering an important simplification occurs. 
As argued in 3], one can then replace the V/j, which are in general off-shell, by their on- 
shell values which only depend on the square of the total cm. energy s = Pq , provided that 
at the same time one uses the physical values of the pion decay constant / — )■ ~ 93 MeV 
and the meson masses (m^r, rriK, Tn^^) ^ (140, 496, 547) MeV that enter into the calculations. 
Then the coupled integral equations for T^^ are replaced by algebraic equations that can be 
solved analytically. In particular, one obtains T/^(s) for KK in both isospin channels by 
first "dressing " the bare KK interaction Vj^i(s) with vrvr or tt^t] polarization loops 1122(5) = 
Il^^[s) or 11,^0^(5) for / = or 1 respectively, and then inserting this dressed interaction into 
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the Lippmann-Schwinger equation for T{^{s) to find the result already quoted in [18 



^ (1 - nin{i)(i - v,iui,) - v,iuiMu{, ^ ^ 

Here n{^(s) = n^^(s) accounts for the common KK polarization loop in both isospin 
channels. 

We characterize the various meson loop diagrams appearing in T/^(s) by the common 
symbol —ill{s) where 

/• rf^/ 1 1 _ 2 

~ J (2vr4) (/2 - ^2) + p^y _ ^2 ' (2) 

that involves the integral over the two meson propagators of masses {nia, mb) in the loop; e is 
a symmetry factor [22|] that takes on the values (1/2, 1) depending on whether rua = mi, refer 
to identical mesons or not. Since this integral diverges for large 4- momenta, we introduce a 
cut-off A in 0(4) momentum space. The on-shell versions of the V/j{s) together with the 
regulated, closed form expressions for n(s) from which IIt,^{s), Uj^Oj^^s) and n^^(s) may be 
constructed, are listed as Eqs. f l20|) . fl2T]) and fl23|) in the Appendix 



The poles of T(i{s) on the appropriate sheet of the cut complex s plane are identified 18 
with the masses and half widths of the scalars /o(600) (or a) and /o(980) for / = 0, and 
ao(980) for 1 = 1. This sheet structure depends in turn on the analytic behaviour of the 
functions n(s) in the denominator of Eq. ([1]) that develop unitarity cuts along the real s-axis, 
starting at the production thresholds of nn, n^rj and KK mesons at s = (ma + mt)^ = 4m^, 
(m7r0 + m^)^ or 4m|- respectively. The prescription for analytically continuing T(^{s) through 
these cuts onto the relevant Riemann sheet in the lower half plane where these poles lie is 
detailed in [isj. 

Using this information together with the closed forms for polarization functions and 
interaction vertices assembled in the Appendix, we reconstruct the /o(980) pole of T'{\(s) as 
a function of the cutoff in the 0(4) regularization scheme. The results are given in Fig. [T]that 
displays the width Tf^ versus the mass M/^, parametrized by A in order to illustrate their 
sensitivity to cutoff. These calculations assume a Breit-Wigner form ~ {PQ — Mff^+iTf^^/2)~^ 
with Pq = ^/s for T{\(s) in the vicinity of the pole. We use these results to fix A by requiring 
that the real part of the pole coincide with the observed /o(980) mass as quoted in either 



the PDG |23| data tables, [(980 ± 10) - i(20 to 50)] MeV, or the Fermilab E791 experiment 



that gives [(975 ± 3) - i(22 ± 2)] MeV. Then 



Mf, - ^Tf, = [(980 ± 10) - i(23 ;^)] MeV for A = 1.35;[j:i^ GeV (3) 

in the first case, while 

- ^r/„ = [(975±3) -i(26T2)] MeV for A = 1.43 ^ 0.05 GeV. (4) 

in the second. These choices for A also reproduce the observed /o(980) half-widths rather 
satisfactorily. Taken together, this suggests that the appropriate values of the cutoff A in 
the 0(4) regularization scheme should lie within a relatively small window around 1.4 GeV. 
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3. Scattering lengths 



The associated KK scattering lengths of good isospin are given in terms of the Tl^{s) by 
the standard expression [22[ 



^KK = - lim, ^ = (5) 



The scattering lengths complex numbers with Im a]^^ < because of the presence 

of the open tttt and vr^r^ decay channels for J = and 1 = 1. One finds, in units of the 
inverse kaon mass mj^^, that 

flO^^ = 3.530 - 2.077i, a]^^ = 1.651 - 1.416i, A = 1.35 GeV (6) 

or 

a^^ = 3.303 - 1.713i, a]^^ = 1.655 - 1.22H, A = 1.43 GeV (7) 

for the two A's found in Eqs. ([3]) and (jlj). 

The values calculated above for the isoscalar scattering length are similar in order of 



magnitude and sign to early direct experimental measurements [25j that yield a^^^ = [(3.13± 
0.30) — (0.67 ± 0.07)z]m^^, and a later analysis [261] of Tin data that gives = (4.36 — 
1.49i)m^^. No similar measurements are available for A model-dependent estimate 



27| extracted from more recent pp — > dK^K^ data [28] yields a]^^ ~ [—(0.05 ± 0.05) — 
;i.59±0.60)]m;^^ These numbers may also be compared with the rough estimate Ur. 



KK 



a]^^ ~ (2.98 — 2.13i)m^^ from the zero-range universal approximation |29(] M — ^ 
2mK — T^K^io.KR)'^ common meson masses and half-widths ~ (980 — 35i) MeV in both 
isospin channels. 

Moving from T-matrices labelled by good isospin to particle-antiparticle labels with the 
aid of Eq. 0181) . one obtains the common strong interaction scattering lengths for the physical 
K^K~ and K^K^ channels as 

aK+K- = ttKORo = -tex + a]^ji) = 2.591 - 1.747i or 2.479 - 1.467i (8) 

for cutoffs (1.35, 1.43) GeV respectively, if the K'^K'^ mass difference A = ttiko —mK± ~ 4 
MeV is ignored. If not, then the K^K~ scattering length becomes 30 



^ / a^+K- fco<^4A- \ ^ 2.688 - 1.896z or 2.567 - 1.566z (9) 



^K 

for the two cutoffs in question. Here ko = y/2mKoA . 



4. Strong interaction energy shifts and decay widths of kaonium 

The level shifts and decay widths for kaonium due to strong interactions have been 
discussed in detail in js], 11] using the vector meson exchange model Lagrangian density 
mentioned in the Introduction. We re-evaluate these quantities using the basic xPT £2 
density given by Eq. ( 1T7|) instead. 
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The unshifted ground state of kaonium lies at Eis = —^a^ji = —6.576 keV where /i = 
\mK± is the reduced mass with = 494 MeV and a ~ 1/137. Contrary to the case of 
pionium, however, where charge exchange tt+tt" — )■ 7r°7r° dominates the decay [l|, K^K~ 
K^K^ is not allowed due the K^K^ mass difference. Thus the principal strong decay modes 
are kaonium — t- tttt + -k^t] that proceed via strange quark annihilation. 

Isospin is broken by the coulomb field as well as the meson mass difference. Since the 
binding energy of kaonium is small relative to that of the strongly interacting KK ground 
state of /o(980) at ~ 10 MeV, the modified energy spectrum can be found by the standard 
procedure [sH of joining the zero momentum s-wave scattering function of mixed isospin 
u{r) = rip{r) — )■ 1— r/op of the pair emerging from the strong interaction zone with scattering 
length Op, onto an exponentially decaying pure coulomb wave at infinity. This is given by an 
incoming Whittaker function [32| PFj^ i/2(2ifcr) at complex momentum k = k\ = —ifiaX with 
ReA > 0. The coulomb parameter index for attractive fields then reads irj = —ifia/k = 1/A. 
Here A = A„, n = 1,2,- ■■ is a set of eigenvalues to be determined by the logarithmic 
matching condition a.t r = d outside the strong interaction zone of extent ~ l/rriK- One 



thus retrieves the Kudryavtsev-Popov eigenvalue equation j33| which we rearrange as 

.A. 1 



a 



pc 



2/i«[^ + lnA„ + ^(l-— )+7], n = l,2,--- (10) 



2 " A, 



Here ip is the standard digamma function 3^, 7 = 0.57721 ■ ■ ■ the Euler constant and a 



'■pc 



1/apc is the physical K^K inverse scattering length in the presence coulomb interactions 
as defined by Bethe [3]| . 

Oipc = otp — 2/ia[ln(2/ia(i) + 7] (11) 

with Up = l/op from Eq. (Q. The role of ape as the relevant experimental observable has 
also been stressed in js^]. We estimate d ~ 2.2m]^^ for the (calculable) Bethe joining radius 



31[ for the coulomb plus strong interaction field; d is thus fixed and not a parameter. Since 
d is much smaller than the Bohr radius 1/fia of kaonium, the small argument approximation 
for Wiri,i/2i2ikr) suffices for deriving Eq. (fTOi) . 

The revised s-wave energy levels and decay widths (A„ is complex because ap and therefore 
ape is complex) of the kaonium atom are given by 

Ex. - -r.„ = ^ = (12) 

In the absence of any strong interactions — )■ so that ape — )■ 00. From Eq. (fTO!) this 
means that A„ is then determined by the poles of ip{^ — j-) as A~^ = n, and one recovers 
the pure coulomb spectrum from Eq. f|T2|) . 

The input for ape on the left hand side of Eq. ffTOj) is obtained from Eqs. ([9]) and ffTTl) as 



a. 



pc 



{2.589 - l.Q5Ai)-^mK or (2.459 - 1.375z)~W (13) 



for A = 1.35 or 1.43 GeV respectively. Eq. ( ITO|) has multiple roots. Solving for the eigenvalue 
Ai corresponding to the Is ground state of kaonium, one finds Ai = 0.9812 + O.OllSi or Ai = 
0.9822 + 0098i for the above two values of the cutoff. The energy shifts AEis = (Ex^ — E^), 
half-widths ^Tu = ^Fai and resulting lifetimes tu = h/Vis (in units 10~^^s) introduced by 



the strong interactions are then [42 
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A^i, - '-Tu = [(246^^1) - ^(304+^^4)] ^ 2.2 t 0.9, A = 1.35;°;}^ GeV 

= [(2331?) - ^(254128)] eV, = 2.6 ^ 0.3, A = 1.43 ^ 0.05 GeV 

(14) 

Notice that the energy shifts are repulsive and of the same order of magnitude as the decay 
widths. 

The uncertainties on A have been also included in these calculations. The larger error bars 
at the lower cutoff clearly overlap with those at the upper cutoff. A conservative estimate 
for the lifetime is thus provided by the spread of values at the lower cutoff of tis ~ 1 to 
3 X 10~^^s, rounded to the nearest integer. 

Fig. |2] gives a parametric plot of Fi^ versus AEi^ for a range of A's. The inset in this 
figure shows an analogous plot of the quantities —Im ape versus Re ape that provide the 
input for generating the main curve. One sees that the two curves behave in a similar 
fashion, and moreover turn at approximately the same value A = 1.18 GeV where Re ape 
reaches its maximum value in the inset of Fig. |5J This behavior is confirmed by the original 
Deser et al. approximate formula j35| that is recovered by expanding the roots of Eq. f lTU]) 
to lowest order in the interaction parameter 2/iaapc about their limiting values = n. 
Then 



AEu - -Tis ^ 2f/a^ape = — apc|V^is(0) |' (15) 

for the Is ground state, where ipisi^) = (/ia)^'^^vr~^/^. In this limit one sees that the two 
parametric curves in Fig. [2] are related by a simple rescaling of their axes. 
The Deser estimate, Eq. (ITSl) gives 

AEr, - ^Fi, ^ [(2481}^) - '-{318+11')] eV, = 2.1 t 0.9, A = 1.35;°;}^ GeV 

^ [(2361^) - ^(2641^0)] eV, = 2.5 ^ 0.3, A = 1.43 ^ 0.05 GeV 

(16) 

for the energy shift and width in the place of Eq. f|T^ . In the present instance this approx- 
imation is seen to be fairly reliable, overestimating the energy shifts and widths by ~ 1% 
and ~ 5% respectively, when compared with the eigenvalue solutions obtained from the 
Kudryavtsev-Popov equation. 



5. Summary and conclusions 

We have studied the energy shift and decay width of the Is ground state of kaonium using 
interactions taken from chiral perturbation theory to construct the strong KK scattering 
amplitudes. These calculations contain a single regulating cutoff A that is fixed by requiring 
that the relevant pole of the KK scattering amplitude in the isoscalar channel reproduce 
the experimental /o(980) mass determinations within their quoted error bars 23, 2J . 
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The calculated decay lifetimes for kaonium then range from ~ 1 to 3 x 10~^^s. Direct 
experimental information on the energy shift and decay width of kaonium is not yet available. 
However, such experiments have already been performed in the case of kaonic hydrogen 
{K~p) by the DEAR collaboration j36|. These data have been analyzed u sing a variety of 
theoretical approaches to extract scattering lengths, see for example [2I, [s^, 38 



A similar experimental program for kaonium would allow one to extract a value for 
the physical K^K~ scattering length by combining Eqs. fll2p and fllOp to determine 
a value for Ai and hence ape- This information would also allow a re-assessment of the 
various theoretical models mentioned in the Introduction for the /o(980) and ao(980) scalar 
mesons by comparing with their predictions for the strong scattering lengths and hence 
Qpc- The additional feature [39^ that kaonium is the only purely mesonic atom with hidden 
strangeness that decays via strangeness annihilation makes such experimental and theoretical 
investigations of particular interest. 

Indirect methods of observation of the formation and decay of kaonium [3I, HqI 41 | may 
also be feasible in the future. 
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6. Appendix 

(i) Lagrangian density C2 and the A-point interaction vertices 

The leading order xPT interaction Lagrangian density for the pseudoscalar meson octet 



used in 18 



IS 



(17) 



where 







7r+ 




* = 


TT 















is a 3 X 3 matrix of the meson fields in SU{3) flavor space constructed from the Gell-Mann 
A matrices that are the generators of this group; r] is identified with r^g. M is the diagonal 
mass matrix M = diag(m^, m^, 2m|- — m^) and / is the pion decay constant. The trace 
runs over the SU{3) flavor space. 

Call iV/j{s) the set of 4-point vertex diagrams generated by iC2- Working in a set of basis 
states of good isospin 



\{KK) 



0,1 



'V2 



(18) 
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and 



|(7r7r)°>= ^[7r+7r- + 7r-7r+ + 7rV], |(A)i>=A (19) 

v3 



the on-shell values j43j for the V/j are 

= — K--f~('^s-^^,nl-\,n^-lml-), V^^'-'f (20) 
Here ^/s is the total coUisional energy in the CM system. 



(ii) Polarization loop integrals 



The expression for the 0(4) regularized integral in Eq. (E]) depends on where s lies relative 
to the cut that starts at the branch point {rria + rribY. ^ > i^a + w^b)^ on the upper lip 
of the cut along the real axis, the integral acquires an imaginary part and one finds 



His) 



(47r)2 

with Lab{s) given by 



ml 



ln(l + 



A2 



mt 



ml 



mt 



1^(1 + ^) - ^'^^(^) 



(21) 



r f ^ 1 1 f ml + ml 

Labis) = -1 - -[— 2 

2 — m^ 



m„ 



mi 



In 



m: 



mt 



tanh ^( 



ah 



+ tanh ^( 



aft 



1 - 



S > ima + mft) 



(m„ - mfc) 2x1/2 / (m^ + m^) 2x1/2 2p, 



(22) 



and Pab is the magnitude of the 3-momentum of either meson in the cm. system. Note 
that all expressions are symmetric under the interchange — )■ m^. Equation fl^T]) with 
rria = m^ and m^ = m^ gives the closed form for n^o^(s). 

For equal masses m^ = m^ = m, n(s) reduces to the simple form 



ms) 



e r , , A' 
l + ln(l + 



(47r)2 

4m2x 1/2 



m^ 



2 2 

m , m , 1 

^(1 + ^)- 



2J/tanh"^ Jf + mJf 



s > 4m 
(23) 



Setting m = m^, or m^ with e = 1/2 or 1 respectively then leads to closed forms for n,r7r(s) 
and IIkx{s). 
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FIG. 1: Mass-width relation for /o(980) parametrized by the 0(4) cutoff A. The rectangle gives 
the limits on Mf^ and Tf^ suggested in the PDG data listings 23|, while the cross indicates the 
Fermilab E791 measurements [2J] and the associated error bars. The KK threshold at 2mK ~ 992 
MeV is indicated by the arrow. The physical parameters used in this evaluation are (m^, m^, /jr) = 
(140, 496, 93) in MeV. 
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FIG. 2: The Is decay width Fi^ versus energy shift AEis for kaonium parametrized by the 0(4) 
cutoff A. The estimated energy shift and corresponding width can he anywhere along the arc of 
the curve intersecting the calculated uncertainties-on-A rectangles centered at A = 1.35 or 1.43 
GeV (large and small rectangles) respectively. These rectangles correspond to fitting either the 
PDG 23|] or the Fermilab E791 data [2j] for the observed /o(980) mass within the experimental 



uncertainties, refer Fig. [TJ The inset shows the —Im ape versus Re Upc parametric curve in units 



m 



^ calculated for the same set of A's. The shape of this curve mimics that of the width-shift 
curve. The latter also turns in the close vicinity of A = 1.18 GeV where the parametric ape curve 



K 



turns as Re ape reaches its maximum value in the inset. 
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